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1.  Introduction 

In  (his  paper,  we  are  concerned  with  the  analysis  of  a  system  operating  in  a  hard-real-time  environ¬ 
ment.  In  such  a  system,  each  job  (or  task)  must  be  completed  within  a  specified  period  of  time  after  a 
request  for  its  execution  arrives.  If  any  job  fails  to  complete  within  its  deadline,  the  entire  system  is  con¬ 
sidered  to  have  failed.  For  a  description  of  such  systems  and  their  analysis  in  a  deterministic  environment 
see  [9j. 

We  consider  here  a  single  server  queueing  system  in  which  job  arrival  stream  is  Poisson  and  the  job 
service  requirement  is  generally  distributed.  Each  job  has  a  deadline  associated  with  it  so  that  if  the 
response  time  of  a  job  exceeds  its  deadline,  we  will  assume  that  the  system  has  failed.  The  system  can 
also  fail  due  to  a  breakdown  experienced  by  the  server.  The  completion  time  (or  the  actual  service  time) 
of  a  job  once  scheduled  will  be  allowed  to  depend,  in  general,  on  the  job  sequence  number.  In  this  way, 
graceful  degradation  of  the  server  can  also  be  taken  into  account.  The  aim  of  this  paper  is  to  derive  an 
expression  for  the  average  number  of  jobs  completed  before  system  failure. 

M/G/l  queueing  system  with  server  breakdown  and  repair  has  been  analyzed  by  Gaver  [4]  while 
and  M/M/n  queueing  system  with  server  breakdown  and  repair  was  analyzed  by  Mitrany  and  Avi-Itzhak 
[11].  Baccelli  and  Trivedi  [3]  analyzed  and  M/G/2  standby  redundant  system  with  breakdown  and  repair. 
These  studies  carried  out  steady  state  analysis.  Approximate  transient  analysis  of  such  a  queueing  system 
has  been  carried  out  by  Meyer  [10]  assuming  no  repair  while  Kulkarni,  Nicola,  Trivedi  and  Smith  [8]  have 
extended  this  analysis  to  allow  for  (possibly  imperfect)  repairs.  The  latter  effort  [7,8],  also  allows  for 
deadline  constraints  to  be  imposed  but  only  in  the  case  that  no  resource  contention  is  permitted.  Analysis 
in  the  present  paper  is  exact  and  allows  for  deadline  constraints,  queueing  and  server  breakdown. 

After  defining  the  basic  model  in  the  next  section,  we  then  derive  the  functional  equations  for  the 
queueing  system  in  section  3.  These  equations  are  specialized  to  the  case  of  an  M/G/l  queue  with 
exponentially  distributed  deadlines  and  independent  geometric  server  failure  process  in  section  4. 
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2.  The  Basie  Model 

Although  various  computational  assumptions  will  be  introduced  later  on,  we  first  state  the  queueing 
problem  in  its  generality. 

Let  {o%  ,n  >  1}  and  {5.  ,n  >  1}  represent  the  respective  completion  time  and  deadline  length  of 
the  n  -th  job.  <r,  represents  the  actual  time  job  n  spends  in  the  server  once  scheduled.  This  may  take 
into  account  a  possible  degradation  of  the  server’s  capacity.  Therefore  the  statistics  of  <rn  will  not  be  star 
tionary  in  general  (for  instance  the  time  to  serve  a  unit-customer  may  increase  with  the  age  of  the  system) 
and  <r%  is  allowed  to  be  infinite  (for  instance  when  the  server  has  failed  before  the  completion  of  the  job). 
Let  t%  ,  n  >  1  represent  the  time  between  the  arrivals  of  customers  n  -1  and  n  .  rn  will  be  assumed  a.s. 


finite. 


We  shall  represent  the  successive  response  times  by  a  sequence  of  defective  random  variables  Rn , 
n  >  0: 

-  R  o  is  a  given  a.s.  finite  random  initial  condition, 

~  ^«+i  represents  the  response  time  of  the  n  +l-st  customer  provided  its  completion  time  is  finite 

and  provided  the  n  -th  customer  experienced  a  finite  response  time  and  met  its  deadline.  Otherwise, 
Rn+i  is  infinite. 

Mathematically,  the  evolution  of  R „  is  described  by  equations  (1)  and  (2): 


{f?.  <oo}f|{<T,+1<oo}p){i?.  <$„} 

{r»+i<  oo}  = 


{/?  1  <  oo}  =  {<7,  <  oo}. 


n  >  1, 


and  in  the  event  {/?,+l<oo}: 

/?,+!  =  [/?.-r,+i]+  +  <th  +,,  n  >0  (2) 

where  [a]+  denotes  max(a,0).  Hence,  on  oo},  RH+l  is  fully  determined  from  (2)  and  the 

knowledge  of  an  initial  value  of  R0.  In  the  complementary  event,  {Rn+1  =  oo},  Rn+l  is  determined  by 


BBS 


We  introduce  now  computational  assumptions  to  be  used  in  the  sequel  of  the  paper: 


-  The  sequences  {r,  ,n  >  1},  {<%  ,n  >  1}  and  {£»  ,n  >  1}  will  be  assumed  independent. 

-  The  arrival  stream  will  be  assumed  Poisson  of  intensity  X. 

-  The  hard  deadlines  will  be  assumed  to  be  i.i.d.  exponentially  distributed  random  variables  with 
parameter  8. 

•  The  <r»'  s  will  be  assumed  to  be  independent  but  not  necessarily  identically  distributed  random 
variables.  We  shall  denote  as  q.  (a )  the  quantity  E  [e  1(0<  <oo)l.  /?«(*)  >  0. 


3.  The  Functional  Equation 

Let  »  GC ,  Re  (e)  >  0.  From  equations  (1)  and  (2)  we  obtain: 

<  o=}l{8.  <  «.)!(<-...  <  co) 
Let 


(3) 


*„(«)  =  £■  1« 1(B.  <0„, |,  »>0.  (4) 

Taking  the  expectation  on  both  sides  of  (3)  and  using  independence  and  distributional  assumptions,  we 
get 


)  =  1.+i(«  JBIe-1*--1-*'1*!,*,  <  Ml(*.  <  „,,)] 

=  <  oo)] 


So  that  4>%  {»  )  satisfies  the  recurrence  equation 


i+i(<)== — — [^»  (*  +5)-*  d,  (X+5)],  n  >  0. 


(5) 


For  z  GO ,  |  z  |  <  1,  let  G(»  ,z )  and  F(»  ,z)  be  defined  by  the  power  series 
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G(a,z)  =  > 

»>i 

F(».*)=  EM*)**  ■ 

•  >o 

Using  the  recurrence  equation  (5),  we  get 

F(s,2Ho(«)  =  T^-E2*+1'lf»+i(a)(X^»(,+5)-a^»(X+5))  •  (6) 

A_4  »  >o 

Then  using  Hadamard’s  theorem  [5]  we  conclude  that  F(«,z)  satisfies  the  functional  integral  equation 
(see  the  Appendix): 

F(s  ,* )  -  Ms  )  +  f  F(»  +6,u  )G  {s  +S±)i*  (7) 

-  ■57—  /  F(X-f5,«)<7(X+tf,-^-)d«]  , 

6i  n  p  » 

where  the  contour  integral  is  taken  on  any  circle  with  center  at  0  and  of  radius  R  where  |  2  \  <  R  <1. 

Before  considering  special  cases,  let  us  note  how  to  get  features  of  practical  interest  from  the 
knowledge  of  the  solution  of  (7).  Let  N*  be  the  stopping  time  defined  by: 

N*  =  inf{n  >  1  |  R.  >6.}.  (8) 

The  very  definition  of  («  )  and  (1)  entail: 

4>. (0)  =  P[F„  <  00]  =  P[ fj (Rk  < 6„ )f)K  <  00))  (9) 

*— 1 

=  a.  P[N*  >  n),  n  >  1  , 


where  a,  ^=P(ffn  <  00). 

So  that  the  mathematical  expectation  of  N*,  the  number  of  customers  served  before  (and  including)  the 
first  hard  failure  is  given  by  (in  case  at  =  a  for  all  n  >  1): 
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kx.- 


E  [N* ]  —  -i-  X  lim^F (0,z )  -  <^o(0)).  (10) 

More  generally,  denoting  by  T  a  circle  of  center  0  and  radius  R  <  1,  one  gets  from  (6)  and  (9) 


(11) 


a2»jr  f.  z' 

4 ■  A  Special  Case:  Systems  with  geometrically  distributed  catastrophic  failures. 

Consider  a  system  in  which  the  server  is  subject  to  catastrophic  non-repairable  failures  (i.e.,  when 
the  failure  occurs,  the  server  stops  functioning  forever).  Assuming  that  the  number  N  of  customers  that 
this  system  can  process  before  the  first  failure  (and  disregarding  hard  deadline  constraints)  is  geometri¬ 
cally  distributed  with  parameter  a  so  that 


P 


P  [N  =  n  ]  =  (l-of)ar*  .  (12) 

Such  a  situation  will  occur  in  case  the  server  can  experience  a  catastrophic  failure  with  rate  \j  so  that 

0° 

the  probability  of  a  successful  job  completion  a  —  /  e  f  dFs(t  )—G  (\{  ), 

o 

where  we  assumed  that  the  service  requirements  of  customers  are  i.i.d.  with  common  distribution  func¬ 
tions  Fg  and  Laplace  Stieltjes  transform  G  (« ).  In  this  case,  we  get  the  following  representation  for 

K  (« )= 


^i(*  )  =  or  G(s) 
and  from  (7)  we  get: 

f(«,*)  =  ^(s)  +  ar-£iil  [\F(s  +  6,z)  -  sF(\  +  S,z)\  . 


(13) 


(14) 


Let 


,(*,*)  =  Xo££lil 


\~9 


b(s,z)  =  M»)~  -  F(\+S, z  ) 


Equation  (14)  can  be  rewritten  as 


(15) 


f.V.'A 


F(*  ,2 )  =  a  [a  ,z  )F(g  + 6,z )+  b(a, 2)  ,  (16) 

so  that  for  any  finite  integer  n  and  for  Re  («  )  >  X 

F(e  ,z )  —  4  {«  ,2 )  +  £  (,no  a  (a  +» 6,2  ))b  (s  +k  6,2 )  (17) 

a 

+  ( Il^a  (<  +i i))F(i  +(*  +l)i,i )  . 

Let  us  show  that  Rn  ( a  ,2 )  the  last  term  in  the  right  hand  side  of  (17),  converges  to  zero  as  n  goes  to 
infinity.  For  large  n  , 


|  a(s  +n  6,2)  | 


^i£_Ll  |  G  (*  +n  6)  | 


so  that  the  convergence  of  this  remainder  is  faster  tb&u  1/n  !  : 

|  /?»(«,*)  |  -/„(*,*)<  K(2)4t- 


Hence  in  the  limit,  (17)  yields  the  expansion: 


F(a  ,2)  —  ^(  TI  a(a+i 6,2 ))b(a+k 6, 2)  , 
k-o  *-° 


where 


II  a  (a  +*6,2 )  ==1  . 

» “0 

Using  our  definitions  (equations  (15))  we  get: 

F(a,2)=  S(Xoz)*(.n  M4  +k  g)  (19) 

k  >0  » —0  a-S  -\6  K  ’ 

-  F  (; \+6,z  )oz  {  £  (Xoz  )*  ( Jl  G  )(,  +k6)}. 

We  reman  with  the  problem  of  determining  the  unknown  function  F  (X+$,z ).  For  this,  let  us  multiply 
both  sides  of  (15)  by  (X-«  )  and  take  the  limit  as  a  approaches  X.  Since  F(a  ,z )  has  to  be  analytic  in  *  at 
X,  the  L.H.S.  has  to  vanish  implying  then  the  relation: 


w*. « 


5]  AiLL-  n  G(X+.*)  •  4(X+(Jt+l)«) 

4-s.n  *!  t~l 

F(X+6,r )  =  X ^ - —g-s - 

E^ff-.n  g(x+,*)(x+m 


where 


One  should  notice  first  that  in  each  term  in  the  R.H.S.  of  (19),  there  are  possibly  singularities  other  than 
«  =  X  situated  inside  the  right  half  plane:  s  =  X  -  *  6  for  those  *  >  1  such  that  X-ti  6  >  0  -  if  any  -  is 
such  a  singularity.  Actually,  we  can  check  directly  that  removing  any  of  these  possible  singularities  (as 
we  did  for  »  =  X)  provides  the  same  expression  for  the  unknown  F(\+6,z )  so  that  the  function 


[AS  )  k  >0  K  \ 


y,^-uk 

x-l  (« )*  ,  w  ,  L  C\  *  >0  * ! 

-  *'*(*)(•+*$) - 7 — TT 

*>o  vlffJl 


( (^)^o(^+(*  +1)^) 


^(X)(X+^) 

*>0  *! 


!/*(,)  =  nG(,+i6)—ii 


where 


» 6+s  -X 


is  an  analytic  function  of  («  ,z )  for  Re(«  )  >  0  and  |  z  |  <1. 

Accordingly,  equation  (11)  provides  the  following  expression  for  the  expected  number  of  customers 
served  before  the  first  failure: 


£[JV*]  = 


E  ~lTui>  (°)Mk  «+*)][  E  TT*'*  W(x+*  6)  ]-[  E  (o)*  6\l  E  tt1'*  (x)*>(x+*+*  5)] 


«/*(X)(X+*5) 


5.  Numerical  Results 


Next  we  give  some  numerical  results  based  on  equation  (22).  In  figure  1,  we  have  plotted  E  [N*  ]  as 
a  function  of  1/6  with  X/  =  0.0001  for  the  case  of  an  exponential  service  time  distribution  with  mean  0.5 
and  the  arrival  rate  X  =  1.  We  also  compare  the  numbers  obtained  by  equation  (22)  with  those  obtained 
by  simulation.  In  figure  2,  we  keep  X j  =  0.0001  and  vary  the  service  time  distribution,  which  is  assumed 
to  be  gamma  distributed  with  mean  0.5.  We  vary  the  shape  parameter  ,a0,  of  the  gamma  distribution 
over  the  values  0.5,  1  and  5.0.  We  have  assumed  in  the  numerical  example  that  #0(s )  ==  1. 

6.  Conclusion 

We  have  studied  an  M/G/l  queueing  system  with  server  breakdown  and  hard  deadlines  on  job 
response  time.  Thus,  a  transient  analysis  of  the  system  is  performed  in  order  to  determine  the  average 
number  of  jobs  completed  before  system  failure.  Extensions  of  the  model  in  the  direction  of  a  more  gen¬ 
eral  “server”  with  multiple  processors,  subject  to  failure/repair  type  of  degradation  in  the  sense  of  [7,8]  is 
needed. 
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Appendix 

Let  G(z)  and  F/i^  be  two  analytic  functions  in  the  domain  |  z  |  <  1  with  expansions 


<H*)=  E  i.2"  . 

»  >0 


r(»)-  s  *.»* , 

n  >0 


Let  T  be  a  circular  contour  of  center  0  and  radius  R  where  |  z  \  <  R  <  1  for  a  given  z  so  that 

I*  I  <1 


For  u  in  the  ring  shaped  domain  |  z  |  <  |  «  |  <  R ,  F(v)  •  G (— )  is  analytical  in  «  and  has  the 


Laurent  expansion 


F(«)  •  G(-)  =  £  E  <t>nlkZ  « 


k  „n-k 


U 


a  >0  *  >0 

Hence,  the  coefficient  of  ,j  G  N  in  this  expansion  is  given  by  the  contour  integral 


.  f(«)G(7) 

E  *•%+/  *'+’  =  27 7  / - TP7 - *  • 

•  >o  £tn  r  « 
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Figure  1.  Average  number  of  jobs  completed  before  system  failure. 


Figure  2.  The  effect  of  varying  the  service  time  distribution 
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